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Abstract 

The tiny neutrino masses measured in the neutrino oscillation experiments can be naturally 
explained by the supersymmetric see-saw mechanism. If the supersymmetry breaking is mediated 
by gravity, the see-saw models may predict observable lepton flavor violating effects. In this work, 
we investigate the lepton flavor violating process /r —> ey in the kind of neutrino mass models based 
on the idea of the “lopsided” form of the charged lepton mass matrix. The constraints set by the 
muon anomalous magnetic moment are taken into account. We find the present models generally 
predict a much larger branching ratio of /z ^ ey than the experimental limit. Conversely, this 
process may give strong constraint on the lepton flavor structure. Following this constraint we 
then find a new kind of the charged lepton mass matrix. The feature of the structure is that both 
the elements between the 2 — 3 and 1 — 3 generations are “lopsided”. This structure produces a 
very small 1 — 3 mixing and a large 1 — 2 mixing in the charged lepton sector, which naturally leads 
to small Br{iJ, —> ey) and the LMA solution for the solar neutrino problem. 
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I. INTRODUCTION 


There are two remarkable features of the neutrino parameters measured in the atmo¬ 
spheric and solar neutrino fluxes experiments, he., the extreme smallness of the neutrino 
masses and the large neutrino mixing angles. The maximal atmospheric neutrino mixing 
and the large mixing angle MSW solution of the solar neutrinos (LMA) are most favored by 
the present experimental data|]^. The two features show that neutrinos have quite different 
properties compared with quarks and charged leptons. It is hoped that understanding the 
neutrino flavor structure may give important feedback on, or even be the key to solve, the 
flavor problem, which includes how to understand the observed quark and charged lepton 
spectrum and mixing. 

Great enthusiasm has been stimulated on studying the neutrino physics and hundreds of 
models have been built and published in the recent yearspj. All the authors tried to reflect 
the two notable features in their neutrino mass models. The smallness of the neutrino masses 
is quite well understood now. See-saw mechanism seems the most natural and economic way 
to present tiny neutrino masses. It is interesting to notice that the study of neutrino physics 
may provide us some information at very high energy scale in the see-saw models. 

As for the large mixing angles, it is fare to say that no satisfying and generally accepted 
explanation is found until now; still, much progress has been achieved in this aspect. As we 
know that the neutrino mixing is actually the mismatch between the bases of the charged 
leptons and the active left-handed (LH) neutrinos, the neutrino models depending on the 
see-saw mechanism are generally divided into two classes according to the origin of the large 
mixing angles coming from, i.e., those from the right-handed (RH) neutrino Majorana mass 
matrix Mr, and transferring to the active LH neutrinos through the see-saw mechanism, 

= , ( 1 ) 

where is the neutrino Dirac mass matrix, and those from the charged lepton mass matrix 
Ml. 

See-saw mechanism is usually achieved in a supersymmetric grand unihed model, which 
has the advantage of gauge couplings unification and avoidance of the standard model hier¬ 
archy problem as well. Since leptons and quarks are placed in the same multiplets in grand 
unihed models and their mass matrices are related together, the question “How to keep small 
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quark mixing while having large neutrino mixing” is raised. The first class of models can 
naturally keep small quark mixing because large neutrino mixing is from the RH Majorana 
neutrino mass matrix TTr, which is unrelated to the quark sector. Several models are built 
based on this idea by different choice of the form of and give satisfied predictions^. As 
for the second class of neutrino models, to avoid large quark mixing, Ml is usually chosen 
very asymmetric, or “lopsided”. Then the large LH lepton mixing is related to the large RH 
quark mixing, which is unobservable. A large number of models have been built based on 
this elegant idea[^ H, 

Neutrino oscillation means the neutrino flavor number is broken, which is conserved in the 
standard model (SM) as an “accidental” global symmetry. Lepton flavor number is another 
conserved quantum number as a global symmetry in the SM. However, it is believed that the 
SM is only a low energy effective theory and it must be extended to a fundamental theory 
at the high energy scale. These global symmetries are then broken in most extensions of the 
SM. Especially the non-zero neutrino masses will also break the lepton flavor symmetry and 
lead to lepton flavor violating (LEV) processes, such as fi —>■ ey, /i — e conversion, fi —> eee, 
T —>• yuy, T —>• ey, Z ^ fir and so on. On one hand, the LEV processes are closely related 
to the neutrino oscillations. On the other hand, they provide different information of the 
lepton flavor structure. Combining study of the two kinds of processes may give us more 
comprehensive understanding of the lepton flavor structure, even leading us to distinguish 
the origin of the large neutrino mixing angles. At least, the LEV processes may give strong 
constraints on the neutrino mass models since they are measured quite accurately and the 
sensitivity will be increased by three or more orders of magnitude in the near future [0, P, 

Recently there appear several works on the topic of lepton flavor violation based on 
supersymmetric see-saw models [|3, |^, 1^ [T^. These works either give predictions of the 
LEV processes based on a specific neutrino mass model or give a general analysis based 
on the neutrino oscillation experimental data and the assumptions of the form of the RH 
neutrino mass matrix. For example. Ref. discusses the branching ratio of —>■ ej 
when the RH neutrinos are degenerate or completely hierarchical so that only the heaviest 
eigenvalue is relevant to the form of Under these assumptions the LEV effects are 
directly related to the neutrino oscillation parameters. Actually many viable models do not 
belong to these kinds. A general discussion of LEV for all neutrino models may be very 
difficult. In our work, we focus on the LEV effects of the kind of neutrino mass models exist 
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in the literature based on the idea of “lopsided” form of the charged lepton mass matrices. 
In this kind of models where the large — v-r mixing in neutrino oscillation is actually 
due to the mixing in the mass matrix of the charged leptons, instead of the neutrinos, it is 
thus natural to expect large LFV effects. In one of our previous works we calculated the 
branching ratios of r —py and Z —>• r/i generally in the “lopsided” models, independent 
of the model details and concluded that the process r —>• /ly might be detected in the next 
generation experiments JT^. We will show in this work that the branching ratio of /i —»■ ey 
is already larger than the present experimental limit in typical models of this kind proposed 
in the literature. We then take i?r(/i —>■ ey) as a constraint which must be satished by the 
models and hnd an interesting structure of the charged lepton mass matrix. The notable 
feature of this model is that in addition to the large 2 — 3 element of the matrix the 1 — 3 
element is of order 1 too. This structure can suppress the branching ratio of p —> ey and 
lead to the LMA solution of the solar neutrino problem at the same time. 

Since there is strong correlation between the muon radiative decay and the muon anoma¬ 
lous magnetic dipole moment {g^ — 2), we furthermore consider the constraints on the SUSY 
parameter space set by the recent BNL E821 experiment ||15||. After the correction to the sign 
of the light-by-light term in the theoretical calculations[|^], there is only a 1.6cr discrepancy 
between the measured value of the muon anomalous magnetic dipole moment and the value 
predicted in the SM[|^. 

The paper is organized as follows. The Sec. II explains how the LFV effects are produced 
in the supersymmetric see-saw models and presents some analytic results of our calculation. 
In Sec III we explain the basic idea of “lopsided” models and show the numerical results 
of i?r(/i —>■ ey) predicted by the models. We then present our new model in Sec IV and 
its predictions of Br{fi —> ey). Finally we give a conclusion of our work and discuss the 
implications of the future experimental results of LFV processes in Sec V. 


II. SOME ANALYTIC RESULTS 

A. Origin of LFV and related formnlas 

In the pure SM the lepton flavor is strictly conserved. If the SM is extended with massive 
and non-degenerate neutrinos, as suggested by the atmospheric and solar neutrino fluxes 
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experiments, the LFV processes may be induced, in analogy to the Kobayashi-Maskawa 
(KM) mechanism. However, such processes are highly suppressed due to the smallness of 
the neutrino masses. The branching ratio is proportional to which is hopeless to 

be observed . 

When supersymmetry enters the theory the scene changes completely. The LFV may also 
be induced through the generation mixing of the soft breaking terms in the lepton sector, 
he., the off-diagonal terms of the slepton mass matrices and the trilinear 

coupling Afj. However, the present experimental bounds on the LFV processes give very 
strong constraints on these off-diagonal terms, with the strongest constraint coming from 
^ ey) (< 1.2 x lO-^i)0. 

A generally adopted way to avoid these dangerous off-diagonal terms is to impose the 
universality constraints on the soft terms at the SUSY breaking scale, such as in the gravity- 
mediated or gauge-mediated p0| SUSY breaking scenarios. The universality assumption 
is certainly the most conservative supposition we can make when we analyze the LFV effects. 
Under the universality assumption off-diagonal terms can be induced through quantum 
effects by the lepton flavor changing operators existing at high energy scale, which are 
necessary to produce the neutrino mixing. 

We calculate these off-diagonal soft terms in the see-saw models, starting from the uni¬ 
versal initial values at the GUT scale, by numerically solving the renormalization group 
equations (RGBs).Then we calculate the branching ratio of /i —>■ ey induced by these terms. 

In a supersymmetric see-saw model, the RH neutrinos are active at the energy scale above 
Mij(We will use the same symbol as the RH Majorana mass matrix Mr to represent the 
energy scale). The superpotential of the lepton sector is then given by 


W = + Yi^H,L,Ej + -M^NiN^ + , 


( 2 ) 


where Vv and Yr are the neutrino and charged lepton Yukawa coupling matrices respectively, 
i, j are the generation indices. In general, Yjv and Yr can not be diagonalized simultaneously 
and lead to LFV interactions. The three matrices Vv, Yr and Mr can be diagonalized by 


Yl 

= uIYrUr , 

( 3 ) 


= vIy^Vr , 

( 4 ) 

Mi 

= X^MrX , 

( 5 ) 
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respectively, where Ul,r , Vr^r and X are unitary matrices. The RH neutrinos can be easily 
integrated out one by one on the bases where Mr is diagonal and then the relation (^) is 
recovered. Dehne 

Vd = UIVl , (6) 

which is analog to the KM matrix Vkm in the quark sector. Then Vr determines the lepton 
flavor mixing. We can see that Vr is determined by the LH mixing of the Yukawa coupling 
matrices Yl and Yat and only exists above the scale Mr. The MNS mixing matrix, which 
describes the low energy neutrino mixing, is defined by 

Vmns = UlU,, (7) 


where t/j, is the unitary matrix diagonalizing the neutrino mass matrix Mj, in Eq. (®- 
If Vmns and Vr are both determined by the neutrino oscillation and LEV experiments 
respectively, it is possible to infer whether the large neutrino mixing is coming from the 
charged lepton sector or not. We can furthermore learn some information about the structure 
of the RH neutrino mass matrix. 


The corresponding soft breaking terms for the lepton sector are 

Csoft = - {mlrE*E, - {mlrN^N, 


+ ( R/ih7ii72 + -BM%N:n; + {A^YlTH.UE^ 




+ {A^Yi,fH 2 L,Nj + h.c. 


( 8 ) 


We assume the universal condition at the GUT scale 



o to 

( 9 ) 

ml 

= rn\ = m\ = ml , 

( 10 ) 

A^ 

= A^ = Ao . 

( 11 ) 


The above LEV effects in the superpotential, determined by Ud, then transfer to the soft 
terms through quantum effects and induce non-diagonal terms from the initial universal form. 
This is clearly shown by the following RGE for m|, which gives the dominant contribution 
to the low energy LEV processes. 



+YLmlYl + ttiI^YlYI + EaE\ + + NaN\ 


( 12 ) 
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with Ea = -Yl, Na = -Yn and gi, Mi being the gauge coupling constants and gaugino 
masses respectively. 

In the basis where Yi and are diagonal, an approximate formula for the off-diagonal 
terms of m| applicable for 5m? -C rUg is given by 

pm|)y S3 + + 

« ■ y4(3 + o")mSlogi^ + ... , (13) 

where in the diagonalized Yukawa matrix Y^ only the (3,3) element Yj^^ is kept under 
the assumption of a hierarchical form of Ym- We will use the GUT-motivated assumption 
Yatj ~ Yt as an order estimate. The ‘a’ is the universal trilinear coupling given by Aq = amQ. 

is the mass of the heaviest RH neutrino. The ellipsis dots '• • •’ in Eq. (|1^) represent 
the corrections due to the RGE running below the scale where the heaviest RH 

neutrino is decoupled. If there is large mixing between the hrst two generations and the 
third generation in V^X, for example, the large 2 — 3 mixing, we will have another two terms 
besides the term in Eq. (p!3D, he., 

■ YH3 + ((1 - |(v3.Y)33n log + |(Vl.Y),y log . 

(14) 

The above term is negligible in the case of V^X 1. 

B. Some analytic formulas 

In this subsection we give the analytic expressions for the branching ratio of the LEV pro¬ 
cess /i — 1 - eq and the muon anomalous magnetic dipole moment within the supersymmetric 
framework. 


The LEV decay of muon occurs through photon-penguin diagrams including sleptons, neu- 
tralinos and charginos, shown in FIG. |l[ The amplitude of the charged lepton radiative 
decay can be written in a general form 

M = emiUj{pj)ia^uq''{A"lPL + A"^PR)ui{pi)e>"{q) , (15) 
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FIG. 1: Feynman diagrams for the process > ey by exchanging charginos and neutralinos 
respectively. 


where Pl,r = |(1 T 75 ) are the chirality projection operators. The j represent initial 
and hnal lepton flavors respectively. The most convenient way to calculate and Aji 
is to pick up the one loop momentum integral contributions which are proportional to 
Uj{pj)PL,RUi{pi)2pi ■ e respectively. The neutralino exchanging contribution is 


aP 

Af 


327r2 ^ ^2 COS ^ 

i'OL 

aP {B^A), 


rriyo 

Aa 




rrii 


(16) 

(17) 


where 


Fi{k) 

F2{k) 


1 — Qk + 3/c^ + 2k^ — log k 

6(1 — kY 

1 — /c^ + 2/c log k 
{1-kf ’ 


(18) 

(19) 


with kaa = Tn?o/irY . A and B are the lepton-slepton-neutralino coupling vertices given by 

Xa ' la 


A 


iaa 


^iaa 




7(iH-3)ci YSa 


ZfiZY + Zp cot Ow) - cot OwYJ - 

^ Mw cos {3 ^ 


2zf^^'^°ZY* + cotew- 


Mw cos (3 


ryiOL ryZa^ 


( 20 ) 

( 21 ) 
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with being the 6x6 slepton mixing matrix and Zjq being the neutralino mixing matrix 
given in Ref. [^, . The corresponding contribution coming from exchanging charginos is 

1 


A? = 


si _ 

327r2^" m? 




ry— ry— * TnjTTlj 

2M2 cos2 13 


m. 




y/2Mw cos P 

92 ^ia* a 1 


^la^2a 


_mj 


rrii 


Fi{ka 


327r2‘"^ m? 


Zt,Zt:F^{Ka) + 




m. 


where 


F,{k) = 

F^ik) = 


V^Mw cos f3 

2 + 3/c — Qk^ + k^ + Qk log k 

6(1 — kY 

3 — 4fc + + 2 log k 


zt: z^: F,{k.,) 


( 22 ) 

(23) 

(24) 

(25) 


(l-fc)3 

with kaa = and Z~ are the mixing matrices of sneutrinos and charginos 

as given in Ref. MM- 


The branching ratio of /x —*• 67 is given by 


a, 


BRii^ 67) = + |/l),T)/r 


dvr 




(26) 


with T^ = 2.996 x IQ-^^GeR 


To identify the dominant contribution one may use the mass insertion approximationlTO 


23| . Under this approximation and when tan/3 is large we have 

'-3 [(mp 


Rr(/i —> 67) 


a 

'Gl 


-^tan2/3 , 
m° 


(27) 


where rris represents general slepton mass and means some constants are omitted on 
the right-hand side. From this expression we can clearly see that the supersymmetric con¬ 
tribution to Br{fi —>■ 67 ) is proportional to tan^ (3 and the hrst two generation slepton mass 
mixing. 

In the computation of the LFV branching ratios we considered the constraints on the 
supersymmetric parameter space set by the BNL E821 experiment of the muon anomalous 
magnetic moment |l^ . We give the related analytic formulas in the following. 

The amplitude for the photon-muon-muon coupling in the limit of the photon momentum 
q tending to zero can be written as 


M' = ieu^ <1 7 " ^ ^ 


(28) 
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with being the muon anomalous magnetic moment. The discrepancy is give by IT? 


6a, = 


-^SM 


= 26(±16) X 10 


-10 


( 29 ) 


The supersymmetric contribution to 6 a, comes from the same photon-penguin diagrams 
producing p —>• ey decay by replacing the hnal state by muon. The analytic expressions for 
6 a, from the neutralino and chargino exchanging contributions are 


H”’ = 


and 

■Sof’ = 


m. 


3271^ cos^ dw rrij 

La 


771 ^ 

y2 7^ ^ia 


+ '^Re{A 

m„. 




(30) 


lOvr^ m? 


{Zt:zta + 


m. 




- 

Xa 


-Re{ZtaZ2a)FA{ka 


(31) 


2M^cos/3' V2Mwcosf3^ 

respectively. 

Comparing the expressions for the amplitude of muon decay and muon anomalous mag¬ 
netic moment we hnd a striking resemblance. Because a mass insertion is necessary to 
give correct fermion chirality in Eq. (p!5D, |y4ij| dominates over the contributions. At large 
tan/? a direct relation between the two quantities is then found under the mass insertion 
approximat ion lO 


Br{^ —»• ey) 


19271 ^ 0 ; 


GfK ' " 


X { 6 a,)‘^ X 


(32) 


where n = A]^/Aj^. It is easy to check that Eq. (|32|) is consistent with Eq. (^7]) and 
K ^ { 6 m‘j^)i 2 /m‘l if we notice that Zt"(m?^)“^Z^" —{ 6 m?^)ij/mg with being the 

common slepton mass. Assuming 6 a, is due to the supersymmetric corrections, Eq. ( |3^ 
may give constraints on the numerical results of the LEV branching ratio if the flavor mixing 
parameter k can be predicted in an explicit model. 


III. p ^ ey IN THE “LOPSIDED” NEUTRINO MASS MODELS 
A. Basic idea of “lopsided” neutrino mass models 

In a SU(5) grand unihed model, the LH charged leptons are in the same multiplets as the 
CP conjugates of the RH down quarks. This feature leads to the fact that the mass matrix 
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for the charged leptons is related closely in form to the transpose of the mass matrix of the 
down quarks. The basic idea of “lopsided” models is that the charged lepton and the down 


quark mass matrices have the approximate forms as[^, 

^0 0 0 \ 

Ml ~ 0 0 CT Md ~ 0 



0 0 e mn 


(33) 




with e -C 1 while a ~ 1 and zeros representing small entries. In the case of the charged 
leptons a controls the mixing of the second and the third families of the LH leptons (Here 


we use the convention that the LH doublet multiplies the Yukawa coupling matrix from the 
left side while the RH singlet from the right side), which contributes to the atmospheric 
neutrino mixing and makes the mixing angle datm large, while e controls the mixing of the 
second and the third families of the RH leptons, which is not observable at low energy. For 
the quarks the reverse is the case: the small 0(e) mixing is in the LH sector, accounting for 
the smallness of Vet, while the large 0(a) mixing is in the RH sector, which is not observable. 

The crucial element in the approach is the “lopsided” form of the mass matrices for the 
charged leptons and the down quarks and that Mi being similar to the transpose oi Md. 
The relation Mi ^ M]y can be achieved in grand unified models based on larger groups 
where SU(5) is a subgroup and plays the critical role to give the relation. This relation 
can also be achieved in models with abelian or non-abelian flavor symmetries, no matter 
they base on grand unification or not. Generally there is no such relationship between Mn 
and the up quark mass matrix Mf/ in SU(5). They are usually assumed to be not lopsided 
and give small mixing in the literature. However, we shall use (Mu)^^ = (Mif)^^ in our 
calculation, which is valid in the simplest SO(IO) model. 

As for the large mixing between the hrst two generations in the neutrino oscillation, which 
is assumed to explain the solar neutrino problem, it is usually given by the LH neutrino mass 
matrix Mi^ in the literature p|, ^ P]. Thus the mixing angle 612 of the LH charged leptons 

is usually small. 

B. Order estimate of p ^ ey in the “lopsided” models 

We hrst give an order estimate of the branching ratio of /x —> ey in the “lopsided” models. 
The quantitative results are given by solving the RGBs numerically and presented in the 
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next subsection. 

The starting point is Eq. In the “lopsided” models there can be no large mixing 

in and the RGBs running below in Eq. (^) is ignored. Under the assumption 

Uv 3 ~ Yt, the ratio between the off-diagonal term ( 5 m|)i 2 and the diagonal term, which is 
approximately mg, is determined by (Ud)i 3 , {^ 0 ) 2 ^ 

Vd is defined in Eq. (^. We estimate Vd through an analysis similar to the analysis in 
Ref. PI in which the order of the MNS matrix element f/e 3 is determined. Writing Ul in 
the form 



1 

0 

0 ' 


Cl3 

0 

Sl3 


C 12 

S 12 

0 ' 

Ul = 

0 

C 23 

■S 23 


0 

1 

0 


-S 12 

C 12 

0 


.0 

—S 23 

C 23 ) 


1 —Sl3 

0 

Cl3 ) 


. 0 

0 



(34) 


where Sij = sin Ctj = cos 9ij and Vl in a similar way with the corresponding angles being 
denoted by 9ij, we get the expressions for (Ud)i 3 and (Vd) 23 , 


(Ud)i 3 — “■Si 2 S 23 Ci 3 — S13C12C23C13 -|- C13C12S13 , ( 35 ) 

(Ud )23 = C12S23C13 + C13S12S13 ~ S13S12C23C13 , ( 36 ) 


with S 23 = sin( 023 —^ 23 )- Under the GUT-based assumption that Yat has a similar hierarchical 
structure with that of up quark, we expect that S 13 is of the same order as the corresponding 
mixing angle, 9l^, in the up quark mixing matrix. If we further assume that no accidental 
cancellation exists between the up and down quark mixing matrices, we then get 9l^ ^ 


Vtd ~ O.OOSp^, where Vtd is the 31 element of the GKM matrix. So, S 13 < 0.008. Similarly, 
we have S 23 ^ Vts = 0.04. One typically finds that S 12 ~ = 0.07 and S 13 ~ 

^ ^12 0, |2^. These two relations hold in most viable models, where m,- and 
get their masses from the 2-3 block of the mass matrix. Gonsidering that S 23 ~ C 23 ~ 0(1) 
in the “lopsided” models then we get 


(Ud)i 3 ~ ~Si 2 S 23 ~ S12S23 ~ 0.05 , ( 37 ) 

(Ud )23 ~ -523 ~ —S23 ~ — 0.71 . ( 38 ) 

The (Vd)i 3 and (Vd )23 are actually determined by Ul solely in the “lopsided” models. This 
conclusion certainly depends on the assumption of the form of Yl and Yat. However, it 
is actually correct in most published “lopsided” models 0, Q. We can actually relax this 
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assumption, unless strong cancellation taking place in Eq. (|35|) , we always get that (Vd)i 3 
is 0(0.05) or larger. 

Mji is determined by using the see-saw relation Eq. (2) conversely 

Mr = (39) 


If we assume that has similar spectrum to the up quarks and Mj, gives the large mixing 
for neutrino oscillation, he., Vl^r ~ Ul,r ~ I and Ui, ~ Vmns, we hnd the scale of Mr is 
actually determined by the lightest neutrino mass, with {Mr)^ ^ 2.5 x lO^^GeE. If Yr is 
lopsided and M^, is approximately diagonal, he., Ur ~ Vmns and Uy ~ /, the scale of Mr 
is then determined by the heaviest neutrino mass, with {Mr)^ 4 x 10^‘^GeV. 

Having known the values of the elements (Vd)i 3 , (^ 0)23 and Mr, we can give the order 
estimation of /i ^ eq by a simple calculation. At the large tan (3 case, the most important 
contribution is coming from the second term in Eq. (^3]). Assuming that all the SUSY 
particles are of the scale irtg we then have 

2 


i?r(p —>• eq) 


^em 5 

^ -— ml 
dvr ^ 


( 

V327r2 


({mluV 

^Xa^taZ 2 aFA{ka) 

\ ml j 

\/2M\y 


C ■ 10-^ 


_7 fWOGeVV /tan/?V 


■ 


m.o 


y V 10 


■ 


tan^/J/r^ (40) 


(41) 


with G being around 1 ~ 10. This estimate means that unless all the SUSY particles tend 
to above ITeU and tan/? is not too large the branching ratio should be above the present 
experimental limit. The reason for such a large branching ratio is clear: the “lopsided” 
structure enhances (¥ 0)23 and (Ud) 13 at the same time, as shown in Eqs. (^) and (pSD. 
The branching ratio is thus approximately proportional to S 23 , which is about 2 orders larger 
if taking S 23 ~ in the symmetric case. Then the branching ratio may be at the 

same order of, or be slightly below, the present experimental limit. 


C. Numerical results 

In this subsection we give the numerical results of the branching ratio of /x —> eq in 
the “lopsided” models. When solving the RGEs we only keep the third generation Yukawa 
coupling eigenvalues Yr and Yatj in Yr and Y^ and ignore the contributions from the hrst 
two generations. We take (¥ 0)23 = —S 23 = and (Ud)i 3 = ■S 12 S 23 = 0.07/v^ as the 
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mo (GeV) 


FIG. 2: 5a^ as a function of mo for tan/3 = 5, mi /2 = 150GeF and for tan/3 = 10, mi /2 = 

150, 250GeF respectively. Aq = niQ and fi > 0 are assumed. For fi < 0, is negative with almost 
the same magnitude as the present corresponding values. The horizontal lines represent the E821 
central value and it TIct and ±1.5(7 bounds. 


typical values in this kind of models. We also show the results when Si 2 is as small as 0.01. 
Below Mr we solve the RGBs of MSSM where is absent. The details about solving the 
RGBS are given in Ref. [Q. We present the dependence of the results on the SUGRA 
parameters: mo, mi/ 2 , Aq, tan/3 and the sign of /i parameter. 

We hrst constrain the SUGRA parameter space by the following conditions: (1) elec- 
troweak spontaneous symmetry breaking is produced by radiative corrections and correct 
vacuum expectation value is given; (2) the LSP of the model be the lightest neutralino with 
its mass limit m^o > 45GeU; (3) the lightest chargino is heavier than 94GeU and (4) the 
lightest charged slepton is heavier than SlGeU. 5a^ gives further constraints on the param¬ 
eter space. It always constrains the p parameter to be positive. We hx Aq = tuq throughout 
the calculation, which does not affect the result much, and take Mr = 4.5 x 10^'^GeU. 

In FIG. D we plot (3a^ as a function of uiq for tan/3 = 5,10 and mi /2 = 150, 250GeU 
respectively. The horizontal lines are the B821 central value and its ±1(7 and ±1.5(7 bounds. 


14 














FIG. 3: Branching ratio of // —> 67 as a function of mo for tan/3 = 5, mi /2 = 150GeF. Aq = mo 
and /U > 0 are assumed. 612 is the mixing angle between the 1-2 generations defined in Eq. (|3^. 
The horizontal dotted line is the present experimental limit for /i —> eq, 1.2 x 10“^^ |0]. 

mi /2 = ISOGel/ is almost the smallest value we can take considering the chargino mass 
experimental limit. We find /i < 0 always predicts a < 0 with almost the same magnitude 
as the corresponding predictions for positive /i. 

The formulas for 5a^ we adopted, Eqs. (POD and (0), are different from those given by 
most authors because we use the full slepton and sneutrino mass matrices including the 
off-diagonal terms. The numerical results, however, do not depend on the mixing angles, 
just as expected. This can be regarded as a check for our program. 


15 







FIG. 4: Branching ratio of /r —> 67 as a function of mo for tan /3 = 10 and mi /2 = 150GeF, 250Gey. 
Other comments are the same as that given in FIG. 

In FIG. ^ > 67 ) is plotted as a function of ttiq for tan (3 = 5, with mi /2 = ISOGeG. 

9 i2 represents the 1-2 mixing angle dehned in Eq. ( |5^ (We omit the bar over 6 thereafter). 
We notice that if 6*12 takes the typical value 0.07 in the “lopsided” models the predicted 
branching ratio is far larger than the present experimental limit in most parameter space. 
For tan /? = 5 we can suppress the branching ratio below the experimental limit when we 
take 6*12 as small as 0 . 01 . 

FIG. I plots Br{fi —> 67 ) as a function of mo for tan/3 = 10 and mi /2 = 150, 25QGeV 
respectively. For the typical value 6*12 = 0.07 the predicted branching ratio is always larger 
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FIG. 5: (5o^ as a function of mo for tan/3 = 20, mi /2 = 250Gey, 350GeF respectively. Other 

comments are the same as that given in FIG. 

than the experimental limit. In case of mi /2 = 150GeI^, even we take Ou = 0.01 the 
predicted branching ratio is still too large for mo ^ 470GeId, where falls within the la 
range. When mi /2 is 250GeV, the branching ratio can be below the experimental limit when 
012 = 0 . 01 . 

We display as a fnnction of mo for tan/3 = 20 in FIG. with mi /2 = 250, 350GeId 
respectively. The valne of tuq in the ±la region can now be as large as abont TOOGeV. 

FIG. H is similar to FIG. ^ except that tan/3 = 20 and mi /2 = 250GeG and 350GeI/ 
now. In case of tan/3 = 20, we must take mi /2 as large as 350GeG so that we can get a 
suppressed —»• ey) when 0 i 2 = 0.01 in the ilex region of 6a^. 

In summary, the branching ratio of /r —> ey is sensitive to the SUSY parameters mo, 
mi /2 and tan/3. For the typical value 0i2 = 0.07, in most parameter space the “lopsided” 
models predict a Br{fi ey) greater than the present experimental limit. Only when 0i2 
is as small as 0.01 can the predicted branching ratio be below the present limit. Following 
this consideration, we build a new kind of models which can produce a near maximal 2 — 3 
mixing, large 1 — 2 mixing in the charged lepton sector and extremely small (Ud)i 3 . 


tan/3=20 
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FIG. 6 : Branching ratio of // —> 67 as a function of mo for tan/3 = 20, mi /2 = 250GeV, 350GeV. 
Other comments are the same as that given in FIG. 

IV. AN INTERESTING FORM OF CHARGED LEPTON MASS MATRIX 


We have shown in the last section that the “lopsided” models usually predict a larger 
branching ratio of /i —> ey than the present experimental limit. The large 2 — 3 generation 
mixing of the charged leptons enhances both the elements (Vd )23 and (Vc))i 3 . Furthermore, 
the scale of Mji is lower than that in the case when the large neutrino mixing comes from 
Mj,. Several authors noticed this fact and pointed out that this might imply a universal 
condition for the soft-SUSY breaking terms at the energy scale not much higher than the 
weak scale BHI , such as in the gauge mediated SUSY breaking (GMSB) models. 
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If we do not give up the attractive supergravity models and at the same time insist that 
the maximal atmospheric neutrino mixing is mainly coming from the charged lepton mixing, 
the sole way to suppress the process /i —>• ey is to suppress (Vd)i 3 . We have examined many 
such models with symmetric elements between the hrst row and the hrst column in the 
charged lepton mass matrix. We found that S 12 is always given by ^Jme/ra^ ~ 0.07 with 
a coefficient of order 1. So the model with such structure will generally give a too large 
branching ratio of /i —> ey, as shown in the last section. However, we find (Vc))i 3 can be 
greatly suppressed in a new form of the charged lepton mass matrix with the elements of 
1 — 3 generations asymmetric. The amazing thing of this kind of models is that a large 
1 — 3 element in can lead to a large 1 — 2 generation mixing, which leads to the LMA 
solution to the solar neutrino problem, and a very small (Vd)i 3 at the same time. As shown 
in Ref. the LMA solution is usually difficult to be constructed and most neutrino mass 
models predict a SMA or VO solution for the solar neutrinos. The reason for this is that 
most models try to produce both the large 1 — 2 mixing, with tan^ Ogoi in the range from 
about 0.2 to 0.8, and the small ratio of mass-squared split ~ 1.4 x 10“^ in 

the neutrino mass matrix M^. This is hard to be achieved and sometimes fine tuning is 
needed. However, if large Ogoi is produced in the charged lepton sector it will be very simple 
to obtain the required neutrino spectrum in M^. 

In the basis where the Dirac neutrino mass matrix is diagonal we give the charged lepton 
mass matrix as 

^ 0 (5 cr' ^ 




-h 0 1 - e 


(42) 


\0 , 1 J 

with 6 = 0.00075, cr' = 0.6 and e = 0.12. e is chosen to fit m^/rrir and 5 to fit rrie/m^. cr' can 
always suppress the value of (Vd)i 3 and at the same time predict a large (Vd)i 2 - We then 
get tan^ 012 = 0.45, sin^ 2023 = 0.997 and (Vd)i 3 = 0.0052. These values are approximately 
the corresponding values in MNS matrix, since if we assume has a similar hierarchical 
structure to that of up quark it will transfer to and give very small mixing in 
If we take a' to fit the value of tan^ 0i2 we have two predictions as (Vd)i 3 = 0.0052 and 
023 ~ 71/4:. We can also take tan^ 0i2 and 623 ~ 7r/4 as two predictions and smallness of 
(fA:))i 3 as the requirement. With a large a', a small (Vd)i 3 and large (Vd)i 2 are produced 
naturally without any fine tuning. 
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FIG. 7: Branching ratio of // —> 67 predicted by our model is plotted as a function of mo for 
tan/3 = 5, mi /2 = 150GeF and tan/3 = 10, mi /2 = 150GeF, 250GeF respectively. Other comments 
are the same as that given in FIG. 

It should be noted that we only want to show an interesting form of the charged lepton 
mass matrix, like Eq. , which can suppress the 1 — 3 mixing and enhance 1 — 2 mixing in 
Vd- We do not intend to build a complete model here and the quark sector is not discussed. 
Even though, it should be noticed that a large a' will not lead to a large mixing in the down 
quarks. 

In FIG. 1^, we show the branching ratio of /i —> eq predicted by the above model with 
tan /3 = 5,10. FIG. ^ shows the predicted branching ratio of the model in case of tan (3 = 20 


20 







mo (GeV) 


FIG. 8 : Branching ratio of // —> 67 predicted by our model is plotted as a function of mo for 
tan/3 = 20, mi /2 = 250GeF,350GeF and for tan/3 = 35, mi /2 = 350GeG respectively. Other 
comments are the same as that given in FIG. 

and tan /3 = 35. In most parameter space in these fignres the decay p —> ey has a branching 
ratio close to the present experimental limit and can be easily discovered by the fntnre 
experiments. 


V. CONCLUSION AND DISCUSSION 

In this work, we calcnlate the branching ratio of p —ey in the nentrino mass models 
with “lopsided” textnre of the charged lepton mass matrix, within the snpersymmetric frame- 
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work. The constraints set on the snpersymmetric parameter space by the mnon anomalons 
magnetic moment are considered. If the charged lepton mass matrix has symmetric ele¬ 
ments between the 1 — 2 and 1 — 3 generations, this kind of models will generally predict a 
mnch larger branching ratio of /i —> 67 than the present experimental limit in most SUSY 
parameter space. 

To accommodate the experimental data, strong constraint should be satisfied by the 
charged lepton mass matrix, he., (Ud)i 3 0.01. We then show a new kind of models with 
asymmetric elements between the 1 — 3 generations, in addition to the asymmetry between 
the 2 — 3 generations in the “lopsided” models. The large 1 — 3 element of the charged lepton 
mass matrix can suppress (Ud)i 3 and at the same time predict a large mixing between the 
1 — 2 generations naturally. This feature is most interesting in constructing the LMA solution 
to the solar neutrino problem, which is usually hard to be achieved in the literature [Q]. 

With the improvement of the sensitivity to the LFV processes by three or more orders 
in the next generation experiments0, it is quite possible that the LFV processes will be 
discovered in the near future, r ^ py is a promising process to determine whether there is 
a large mixing between the 2 — 3 generations in the charged lepton sector, such as in the 
“lopsided” models[Q. If both r ^ /ly and p —*• ey are discovered in the near future, a kind 
of lepton structure as given by the present work should be very attractive. If only r —/ly 
is found, this will imply a large mixing between the 2 — 3 generations in the charged lepton 
sector. However, very special structure should be designed to suppress the p —>• ey process 
further. On the contrary, if r —py is not found while /i —>• ey is found, this should be easy 
to understand and in most models without a lopsided structure this case is predicted. The 
last possibility is that no LFV is found at all. In this case the most natural explanation is 
that the soft-SUSY breaking parameters are universal at the energy scale not much higher 
than the weak scale, such as in the GMSB mechanism. 

If LFV is found in the near future, the particles inducing LFV, such as the SUSY particles, 
may be produced directly on the high energy colliders. The measurements on high energy 
colliders will help us to infer the flavor structure when combining with the LFV experimental 
results. At the same time, interplay with the experiments of neutrino oscillation is important. 
If LFV is found and at the same time Uez of the MNS matrix is measured, much more 
constraint is obtained in building the neutrino mass models. In any case, the measurement 
of lepton flavor violation is quite important to infer the possible lepton flavor structure or 
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other possible new physics, no matter the result is positive or null. 
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